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ABSTRACT 

Context. Herbig Ae/Be stars (HAeBe) are pre-main sequence objects in the mass range 2 Mq < M* < 8 M©. Their 

X-ray properties are uncertain and, as yet, unexplained. 

Aims. We want to elucidate the X-ray generating mechanism in HAeBes. 

Methods. We present a XMM-Newton observation of the HAeBe HD 163296. We analyse the light curve, the broad 
band and the grating spectra, fit emission measures and abundances and apply models for accretion and wind shocks. 
Results. We find three temperature components ranging from 0.2 keV to 2.7 keV. The O VII He-like triplet indicates a 
X-ray formation region in a low density environment with a weak UV photon field, i. e. above the stellar surface. This 
makes an origin in an accretion shock unlikely, instead we suggest a shock at the base of the jet for the soft component 
and a coronal origin for the hot component. A mass outflow of Af s hock ~ 10~ 10 Mq yr" 1 is sufficient to power the soft 
X-rays. 

Conclusions. HD 163296 is thought to be single, so this data represent genuine HAeBe X-ray emission. HD 163296 
might be prototypical for its class. 

Key words, stars: formation - stars: individual: HD 163296 - X-rays: stars 



1. Introduction 

The phase between protostar and main-sequence object is a 
key stage for planet formation and, furthermore, the prop- 
erties of a star are determined precisely during that phase 
for the rest of its life. It is therefore important to study ob- 
jects of this age in order to understand the origin of plan- 
etary systems. Herbig Ae/Be stars (HAeBes) are thought 
to be the predecessors of main-sequence (MS) stars in the 
mass range 2 M Q < M* < 8 M Q , although the empir- 
ically denned object class of HAeBes may in fact encom- 
pass a wider range of evolutionary stages such as stars with 
magnetospheric accretion (Ae) and boundary layer objects 
(Be). HAeBes are certainly young, they do have clear sig- 
natures of surrounding disks and thus must be considered 
to be the more massive brothers of the better studied clas- 
sical T Tauri stars (CTTS), which are young (< 10 Myr), 
low mass (M» < 3M Q ), pre-main sequence stars exhibiting 
strong Ha emission. CTTS are also surrounded by disks, 
as evidenced by a strong infrared excess, and are actively 
accreting material. 

Strong X-ray emission is a characteristic property of 
most young stellar objects. A review of the observa- 
tional situation before X MM-Newton and Chandra fo r 
low mass stars is given by iFeigelson fc Montmerlg ([1999). 
Zinnccke r fe Preibischl (|1994f ) were the first to carry out 



* Based on observations obtained with XMM-Newton, an ESA 
science mission with instruments and contributions directly 
funded by ESA Member States and NASA. 
** Fig. [3] is available in electronic form via 



http:/ /www. edpsciences.org 



a systematic X-ray survey of HAeBes. Restricting atten- 
tion to those stars located within 500 pc, 70 % of all 
HAeBes were found to produce detectable X-ray emission 
with X-ray luminosities ranging between a few 10 29 erg s _1 
and 10 31 erg s , thus exhibiting a much larger detection 
rate th an found for field A - type stars (10 -15 %, accord- 
i ng to ISchroder fc Schmittl (|2007f )). Also, iGiardino et all 
(2004) report the detection of a large X-ray flare (with 
XMM-Newton) from the Herbig Ae star V892 Tau, and 
present strong arguments that the flare in fact originated 
from the Ae star in the V892 Tau system; note, however, 
that ISmith et all (|2005h show V892 Tau to be a close bi- 
nary. Until today it is unclear whether the observed X-ray 
emission originates from the HAeBes themselves or from 
unresolved companions, which by ne cessity would have to 
be low-mass, young, and active stars. ISkinner et alj (2004) 
studied a sample of ten close HAeBes and showed that 
the X-ray emission probably originates from a magnetically 
confined plasma, although it remains unclear if this plasma 
is associated with a companion. In a sample of 17 HAeBes 
IStelzer et al.1 (l2006ri find X -ray emission in about 80 % of 
the objects, more than can be reasonably expected from 
late-type companions, although it is possible that their 
sample i s biased towards kn own X-ray sources. In a follow- 
up paper Stelzer et ail (|2008h detect every object, but again 
many of them show spectral characteristics compatible with 
low-mass companions. 

X-ray astronomy has been revolutionised by the avail- 
ability of high resolution grating spectroscopy in the current 
generation of X-ray satellites. Grating spe ctra of HAeBes 
published so far are taken from AB Aur (|Telleschi et alj 
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2007T) and the spectroscopic binary HD 104237 (|Testa et all 
2008) , where the main component is a Herbig Ae star with 



Table 1. Observing log. 



a CTTS companion of spectral type K3. More observations 
are available from CTTS; their high resolution X-ray spec- 
tra typically show a soft component (which may be hidden 
by large absorption columns) and unusually low f/i-ratios 
in the He-like triplets of Ne and O. For CTTS the most 
promising explanation for these phenomena is a magnet- 
ically funnelled infall model. The inner disk is truncated 
at the corotation radius and ionised material is loaded 
onto the field, flowing along the magnetic field lines and 
hitting the stel lar surface close to free-fall velocities (e.g. 
IShu et aTlll994l ) . In the photosphere an accretion shock de- 
velops, which produces the soft X-ray component. Because 
the plasma has relatively high densities, the f/i-ratios in 
the Ne IX and O VII triplets are small. This model has been 
successfully applied in explaining the X-ray emission from 
the CTTS TW H ya in terms of an acc retion shock model 
plus a hot corona (jGiinther et al.ll2007l ). 

A question unanswered to date is whether the accretion 
disks surrounding HAeBes lead to similar phenomena as 
those surrounding CTTS or if other modes of X-ray gen- 
eration operate in the more massive stars. This could be, 
e.g., magnetically confined winds coll iding in the equato- 
rial pl ane as suggested for IQ Aur bv lBabel fc M ontmerlc 
( 1997f , or i nternal shocks in u n stable winds as in the CTTS 
DG Tau (iGiidel et all 120051 l2008t ISchneider fc Schmittl 
l2008t iGunther et al.ll2008Th 

All these models require magnetic fields, which are 
only weak in evolved A and B stars, because of the ab- 
sence of an outer convection zone required for a solar- 
type a — il— dynamo. Nevertheless the HAeBes could have 
magnetic fields, produced by the compression of a primor- 
dial field of the proto-stellar cloud. It has been suggested 
that they are the p rogenitors of the magnetic Ap/Bp stars 
(|Wade et al.ll2005l h which comprise about 5% of the total 
A star population. 

In order to enlighten the origin of the X-rays from 
HAeBes we performed X-ray observations of HD 163296 
with XMM-Newton, whose stellar properties we explain in 
Sect. [21 before presenting the observations in Sect. [H We 
show the results of our analysis in Sect. 0] and discuss their 
implications in Sect. [5] Sect. Ogives our conclusions. 

2. Stellar properties 

HD 163296 is an isolated HAeBe far from any na- 
tal molecular cloud. With its distance of 1221^ pc 
(|van den Ancker et al.1 11998[ ) it is one of the closest ob- 
jects of its kind, making it an ideal target for studies of 
young A stars. The presence of cold dust around HD 163296 
has lo ng been known from spatially unresolved observa- 
tions (|Manningsl Il994h . and a few years ago the disk 
of HD 163296 has been coronographically imaged with 
HST /STIS showing indi cations for a planetary body in the 



disk (|Gradv et al.l 120001 ). Furthermore, these observations 
allow to trace an outflow in the Lya line fr om 7.3 AU to 
725 A U with a velocity of 350 km s" 1 (|Devine et al.l 
2000). This jet, called HH 409, contains knots up to about 
3000 AU from the c entral star, with an a symmetry between 
jet and counterjet (|Wassell et al.l [20061. The surrounding 
disk has also been imaged in the millimetre range indicat- 
ing a strongly evolved disk, wher e larger bodies a lready 
influence the evolution of the gas (|Isella et al.ll2007l ). The 



Observatory 


ObsID 


Obs. date 


Exp. time 


Chandra 


3733 


2003-08-10 


20 ks 


XMM-Newton 


0144271401 


2003-10-11 


7 ks 


XMM-Newton 


0502370201 


2007-09-23 


106 ks 


XMM-Newton 


0502370301 


2007-09-24 


22 ks 



spectrum of HD 163296 displays a strong infrared excess, 
which can be attributed to the inner disk regions. Apart 
from a single outburst the optical light curve has been re - 
markably constant over the last 25 years (ISitko et al.l l2008). 

The star itself is of spectral type Al with an effective 
temperature T c g = 9300 K, a radius R = 2.1 Rq and a 
mass M = 2.3 Mq. According to evolutionary models this 
places it at an age of about 4 Myr. With an Ay = 0.25 
the s t ar is only mi l dly ab sorbed (|van den Ancker et al.l 
I1998D . iHubrig et al.l (|2006[ ) performed spectropolarimet- 
ric observations of circularly polarised light at the VLT 
for some HAeBes, but were unable to detect a signifi- 
cant longitudinal magnetic field (Bz = —57 ± 33 G) on 
HD 163296; they remark, however, that a monitoring cam- 
pain , which could detec t accretion funnels, is still miss- 
ing. iDeleuil et al.l (|2005f ) present Far ultraviolet spectro- 
scopic explorer (FUSE) and HST/STIS observations of 
HD 163296 showing chromospheric signatures and over- 
ionised species. These can be explained by either an ac- 
cretion shock or in a magnetically confined wind model. 
Hydrogen Lya and C in lines show outflow signatures. 
ISwartz et all (|2005t ) looked for signatures of a companion, 
but HST/STIS imagery excludes binarity down to 0'.'05 sep- 
aration. Their long-slit STIS spectra also show only the 
emission of an early- type unresolved point source. 

HD 163296 is a ROSAT X-ray source, detected in the 
ROSAT all-sky survey (RASS) data with an elongated 
source near the position of HD 163296 and resolved in 
three sources with a n unpublished ROSAT /HRI pointing. 
ISwartz et al.l (|2005f) confirm the detection of X-ray emis- 
sion from HD 163296 using Chandra ACIS-S in imaging 
mode. They find an extremely soft spectrum which can be 
well fitted with a 1-temperature model (kT k, 0.5 keV, 
where k is Boltzmann's constant and T is the temperature). 
Additionally there are five X-ray photons at the position of 
the knot H in the jet. This is a significant source detection 
at the 95 % confidence level, with a luminosity more than 
two orders of magnitude below the central component. 



3. Observations and data reduction 

We observed HD 163296 for 130 ks with XMM-Newton with 
the RGS as prime instrument applying the medium filter 
to block out the bright optical radiation. The observation 
is split in two exposures in consecutive orbits. Additionally 
we retrieved archival observations from XMM-Newton and 
Chandra. Table [1] summarises the observation information. 
HD 163296 is optically bright, just below the formal bright- 
ness limits of the UVW2 and UVM2 filters, but all expo- 
sures of the optical monitor are overexposed and unusable 
for analysis in the inner region. All data was reduced us- 
ing standard XMM-Newton Science Analysis System (SAS) 
software, version 7.1 or the Chandra Interactive Analysis 
of Observations (CIAO) software, version 4.0, in the case 
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of the Chandra observation. Because the exposure is par- 
tially contaminated by energetic proton events we applied 
a time filtering, where a good time interval is defined in 
the usual way for the EPIC/PN camera as < 1 cts s _1 . 
For the EPIC/MOS detectors we lowered the cut values to 
0.25 cts s _1 in order to suppress the high energy noise in 
the spectra. CCD-spectra and count rates were extracted 
from a circular region within the 15" around the target for 
the MOS. We obtained RGS spectra and found a strong 
contamination in the RGS1. We therefore reduced the ex- 
traction region for the RGS1 to the central 66% of the point 
spread function (PSF), keeping 90% of the PSF extraction 
region for the RGS2. This step helps to suppress the back- 
ground. To increase the signal we merged the two exposures 
for each RGS using the SAS task rgscombin e. Spectral fit - 
ting was carried out using XSPEC VI 1.3 (j Arnaudl fl996l) . 
and individual lin e fluxes were measured using the CORA 
line fitting tool (jNess &: Wichmarml l 2002V Because the line 
widths are dominated by instrumental broadening we keep 
them fixed at AA = 0.06 A. 

For comparison purposes we also analysed the XMM- 
Newton data on AB Aur (ObsID 101440801) wh i ch wa s 
already presented in great detail bv lTelleschi et al.l (|2007ft . 
We extracted the central 80% of the PSF in order to min- 
imise the contribution from the close source SU Aur and 
measured the line fluxes in the He-like triplets of neon and 
oxygen with the same method used for HD 163296. 

4. Results 

We first present the variability observed in HD 163296, be- 
fore we show fits to the emission measure distribution. Then 
we give individual line fluxes and take a closer look at the 
density and UV-field sensitive He-like triplets. We compare 
the O vni to O vn line ratio with the value found for MS 
stars. 

4.1. Variability 

Figure [1] shows the light curve and hardness ratio of the 
observations taken in 2007. They are binned to 2 ks. 
We define the hardness ratios as a count ratio of (hard- 
soft) /(hard+soft), where the hard band is taken in the 
range 0.8-2.0 keV and the soft band as 0.2-0.8 keV. 
HD 163296 exhibits only modest variability over most of 
the observations with a sudden increase of the count rate 
by 30% about 30 ks in the observation. No correlation be- 
tween the hardness and the luminosity is visible. In the 
beginning low count rates coincide with a hard spectrum, 
but this pattern does not persist. 

To analyse the long-term variability we compare the 
count rates in the EPIC/PN detectors of all XMM-Newton 
observations. In 2003 the average count rate was 0.166 ± 
0.007 cts s" 1 and in 2007 0.202 ± 0.002 cts s" 1 and 
0.218 ± 0.004 cts s — 1 . The errors are statistical and do 
not represent temporal variation within each exposure. 
The exposure in 2003 is too short to extract a meaning- 
ful lightcurvc. We obtained a i?0<S^4T/HRI count rate of 
0.012 cts s _1 for a pointing on HD 163296 observed in 
1995 from the HRI catalogue and converted it to an energy 
flux in the band 0.3-2.0 keV using WebPIMMS assuming 
that a spectral model fitted to the XMM-Newton data is 
applicable. Fitting the same model to the Chandra CCD 
spectrum we also integrated the model flux in the same 
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Fig. 1. Lightcurve and hardness (soft band: 0.2-0.8 keV, 
hard band: 0.8-2.0 keV) during the observations in 2007. 

band as the other observations. The RO SAT /HRI lumi- 
nosity is about 4.3 x 10 29 erg s _1 , in 2003 Chandra mea- 
sured 3.8 x 10 29 erg s"\ with 4.1 x 10 29 erg s" 1 the XMM- 
Newton observation taken in the same year shows a flux 
very similar to the value from 1995. Our new observations 
are slightly brighter, with average fluxes of 4.5 x 10 29 erg s" 1 
and 4.8 x 10 29 erg s _1 respectively. The uncertainty of the 
luminosities is dominated by the intrinsic variation, not 
by count statistics. We find short-term variations of 30% 
within a single exposure, and the separate observation dif- 
fer on a similar scale. Thus, we conclude that the total 
luminosity has not changed much over the past decade in 
accordance with the optical behaviour (Sitk o et"al] [2008). 

4.2. Emission measure and abundance 

Because of the absence of significant variability during our 
observation in 2007 we jointly fit all available spectra from 
EPIC/MOS, keeping them separately, and the two RGS 
spectra, with the two exposures merged for each detec- 
tor, again keeping the two detectors separately. The spec- 
trum can be satisfactorily described by three thermal com- 
ponents with variable abundances. We use VAPEC mod- 
els and give abundances relative to the solar values from 
iGrevesse fc Sauval (|1998| ). Our best fit model is shown in 
Fig. [2 where it is plotted as line on the MOS1 data. The 
data is binned to contain at least 15 counts per bin. We 
show the data only up to 3 keV because the source flux 
drops quickly and the high-energy spectrum is dominated 
by residual noise in the data. The spectrum is nearly flat on 
the low energy side, indicating a weak absorption. Several 
emission peaks are visible, most notably the Lya lines of 
Mgxn at 1.35 keV and Sixiv at 1.85 keV. Line contribu- 
tions identified in the high resolution RGS spectra are listed 
m table H 

In tables [2] and [3] we give the best fit parameters for our 
model and the associated errors (90% confidence interval). 
The reduced x 2 value for our model is only 1.1, but fewer 
components cannot reproduce the data: A fit with two or 
only one component fails to produce the correct slope below 
1 keV. This can only be compensated by extra emission 
from unresolved carbon lines, which in turn requires a very 
low contribution from nitrogen around 0.5 keV. In total 
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Fig. 2. Low resolution spectrum (MOS1) of HD 163296 and 
our best-fit model. 

Table 2. Best-fit model parameters (90% confidence inter- 
val). 



component 


soft 


medium 


hard 


kT [keV] 

EM [10 52 cm" 3 ] 


u - zl -o.oi 

9 q+0.8 

z '°-0.4 


u - J1 -0.03 
1 9+O.5 
1 ' z -0.6 


9 7+1.6 
z - ' -0.8 
U.O_ 2 


7V H [10 2U cm^] 




7 +4 

' -4 





Table 3. Abundance of elements and first ionisation po- 
tentials (F1P); errors show 90% confidence intervals. 



Element 


abundance 


FIP [eVJ 


C 


q 7+3.0 

°-'-1.5 


11.3 


O 




13.6 


Ne 


1 9+0.5 
± - z -0.6 


21.6 


Mg 


9 q + 1.0 
z -°-0.8 


7.6 


Si 


9 q+1.2 
^■ 8 -0.8 


8.1 


Fe 




7.9 



these models end up with a carbon abundance an order of 
magnitude above solar and virtually no nitrogen. We deem 
this an unphysical scenario and conclude that at least three 
temperature components are present in HD 163296. Within 
the error the fit using three components does not change 
if nitrogen is taken as a free parameter, but the error on 
its abundance is so high that we prefer to keep it fixed at 
the solar value. As a cross-check we test the model with the 
data of the EP1C/PN detector, which has a larger effective 
area, but lower energy resolution; for this case we find the 
reduced x 2 cd = 1.1. 

Table [3] shows the abundances as determined from the 
combined fit. Elements with a low first ionisation potential 
(FIP) are enhanced. A similar pattern is typical of inactive 
stars and its signatures c an be found in many X-ray data 
sets (|Robrade et al.ll2008| ). 

We measure the absorbing column density towards the 
source to Ah = 7l 4 x 10 20 cm~ 3 . Assuming a standard 
gas-to-dust ratio the optical reddening and the X-ray ab- 
sorption should be related through the formula Ah = 



A v -2 x 10 21 cm" 2 (|Savaee fe Mathislll979l see lVuong et~aT1 
(2003) for a compilation of other conversion factors in the 
literature, all roughly consiste nt with this value). The op- 
tical reddening of Ay — 0.25 (|van den Ancker et al.l ll9981 
is fully consistent with the X-ray value. 

4.3. Line fluxes 

The spectrum is dominated by emission lines, it is shown 
binned to a minimum of five counts per bin in Fig. [3] (online 
only). Applying the global model to the RGS data alone 
gives Xrcd ~ significantly more than in the joint fit. 
This is not surprising since our model of only three temper- 
ature components is a simplification of the real temperature 
distribution. Lines which are very sensitive to small temper- 
ature differences cannot match the model precisely. The fit 
can be improved to x^cd = 1-1 by increasing the abundance 
of carbon by a factor of two and iron by a factor of four. 
As explained above we used a smaller extraction region for 
the RGS1. A comparison between spectra extracted using 
the standard extraction region and our reduced extraction 
region shows that our choice of extraction area reduces the 
noise level. The line fluxes measured with both methods 
are compatible within the errors. We list the line fluxes in 
table HI Due to zero-point offsets the measured wavelength 
and its theoretical value can differ; in the table we give 
the theoretical wavelength of the identified lines together 
with the fitted values. To fit multiplets we keep the wave- 
length difference between the components constant. Only 
three lines are measured in both RGS detectors. In this 
case we use the error- weighted mean of the RGS1 and the 
RGS2 flux. The photon fluxes of Ovill Lya are fully con- 
sistent between the RGS1 and the RGS2 within the errors; 
for Fexvn 15. 01 A and Fexvn 16.78A the derived fluxes 
differ by more than the formal lcr-error due to the uncer- 
tain determination of the background level in particular in 
the presence of the very close Fexvn 15.26A line. 

We significantly detect lines of neon, oxygen and car- 
bon in the He-like or H-like ionisation stages and, with a 
significance just above 2cr, nitrogen. We also find a number 
of Fexvn lines, which have a peak formation temperature 
around 8 x 10 6 K. From the fitted line fluxes we calculate 
the t otal line intensity, using the absorption c ross section 
from iBalucinska-Church fe McCammonT (|1992T ). The rela- 
tive errors on the line intensity are larger than those on the 
line fluxes, because here the uncertainty in the fitted Ah 
and the distance contributes to the total error budget. 

4.4. He-like triplets 

Of special interest are the line fluxes of the He-like triplets 
in O vn and Ne ix. These triplets consist of a recombina- 
tion (r), an intercombination (i) and a forb idden (f) line 
(|Gabriel fe Jordanlll969t|Porquet et al.ll2001[ ). So-called R- 
and G-ratios are defined as R = f /i and G = (/ + i)/r re- 
spectively; for high electron densities or strong UV photon 
fields the i?-ratio drops below its low-density limit, because 
electrons can be excited from the upper level of the forbid- 
den to the intercombination line collisionally or radiatively. 
Figured] shows the O vn He-like triplet in the RGS1 detec- 
tor. The i?-ratio is 4.6, which is above the low-d ensity ratio 
of 3.4 as obtained from t he CHIANTI databas e dDere et al. 
llIM lLandi et al.ll200ll or 3.95 from APEC dSmith et al. 
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the line width, for the fit we fix the difference in A 



Line ID 


A (theory) 


RGS1 


RGS2 




flux 


unabs. intensity 




r X l 

[A] 


A (fit) [A] 


[counts] 


A (fit) [A] 


[counts] 


[10 b cts cm 2 s x ] 


r-i *a97 — 1 i 

[10 erg s L ] 


Ne x Ly a 


12.14 


n.a. 


n.a. 


12.14 ± 0.01 


38 ± 9 


5.1 ± 1.2 


17.7 ± 6.2 


Neix r a 


13.46 


n.a. 


n.a. 


13.45 ± 0.01 


32 ± 9 


4.4 ± 1.2 


14.5 ± 5.7 


Neix i a 


13.56 


n.a. 


n.a. 


13.55 ± 0.01 


21 ± 9 


2.9 ± 1.2 


9.6 ± 4.8 


Neix f 


13.70 


n.a. 


n.a. 


13.69 ± 0.01 


22 ± 8 


3.0 ± 1.1 


9.9 ± 4.6 


O viii Ly a 


18. y ( 


18.98 ± 0.01 


74 ± 10 


i o no j_ n m 
18.98 ± 0.U1 


I 1T_L 19 

II / ±13 


18. ( ± 1.2. 


c a a _i_ on a 
54. D ± 20.4 


O VII r 


91 (\ 
Z1.0 


Zl.o ( ± U.U1 


4z. (±1.0 


n.a. 


n.a. 


1/1 (\ -t- 9 1 
14. D It Z. t 


A K 4- 91 ^ 


O vii i 


21.8 


21 77 ± 01 


7.8 ± 4.4 






3.2 ± 1.8 


1 o n + 7 9 


Ovu f 


22.1 


22.07 ±0.01 


46.2 ± 8.1 


n.a. 


n.a. 


14.6 ±2.6 


46.0 ±22.6 


Nvn Lya 


24.78 


n.a. 


n.a. 6 


24.78 ±0.01 


14 ±6 


3.3 ± 1.4 


8.0 ±4.6 


C VI Ly a 


33.70 


n.a. 


n.a. 


33.70 ±0.01 


50 ±9 


20 ±3.6 


73.0 ± 56.6 


Fe XVII 


15.01 


15.01 ±0.01 


50 ±9 


15.03 ±0.01 


101 ± 13 


11.5 ±1.7 


35.5 ± 11.6 


Fe xvn 


15.26 


n.a. 


n.a. 6 


15.28 ±0.01 


68 ± 11 


8.7 ± 1.4 


26.8 ± 9.0 


Fe xvn 


16.78 


16.79 ±0.01 


26 ±7 


16.78 ±0.01 


59 ±11 


6.6 ± 1.3 


20.0 ± 7.6 


Fe xvn 


17.05 d 


n.a. 


n.a. 


17.05 ±0.01 c 


53 ± 14 


7.4 ± 1.9 


22.4 ±9.3 


Fe xvn 


17.10 d 


n.a. 


n.a. 


17.10 ±0.01 c 


59 ± 14 


8.2 ± 1.9 


24.9 ± 10.0 



a line blended (see text) 

6 region contains empty bins 

c lines are not resolved, but two components are necessary to reproduce 




21.4 21.6 21.8 22 22.2 

A / A 

Fig. 4. O vn He-like triplet with our best fit to recombina- 
tion, intercombination and forbidden line. 

I2001I L The error on this ratio is dominated by the statisti- 
cal error on the weak i line and is significantly asymmetric. 
We run a Monte-Carlo simulation to obtain lower limits on 
the i?-ratio, which are described in detail in appendix [XI 
We determine the background from the line-free regions in 
the range 20-25 A to 30 counts A -1 . At a 90% confidence 
level the i?-ratio is higher than 2.6, at a 99% confidence 
level the lower boundary is 1.7, it is thus fully compati- 
ble with the low-density limit, but not with a high-density 
case. The G-ratio is 1.21q4 (90% confidence), where the 
errors are based on the same method. According to the 
CHIANTI database this diagnoses a plasma temperature 
between 0.02 keV and 0.2 keV, a range below the coolest 
component fitted in the global emission measure analysis, 
but in nature the plasma is not divided between three com- 
ponents of fixed temperature, instead each component has 
to be interpreted as a representative for a range of temper- 
atures. Ovu forms only on the cool end of that distribu- 
tion, above 0.20 keV O viii is the dominant ionisation stage. 



13.2 13.4 13.6 13.8 

A / A 

Fig. 5. Neix He-like triplet with our best fit to r, i and f 
line. The i lines is blended with strong iron emission. 

Therefore it is not surprising to find a lower temperature 
from the O VII G-r atio. A small G-ratio cou ld also be due 
to photoexcitation (jPorquet fc Dubaul l2000). 

The interpretation of the Neix triplet is more difficult 
as the i line can be strongly blended by iron lines, predom- 
inantly Fexix and Fexx. The signal is not suffi cient for a 
detailed deblending as in lNess fe Jordan! (|2008f ). However, 
there are two indications, that the counts in the i line are 
dominated by iron contamination: First, we fitted the He- 
like Ne triplet keeping the relative wavelength constant, 
the results of this step are given in table |4] and shown in 
Fig. [5] Then we kept the wavelength for the r and f line 
where the blending is less severe and set it free for the i 
line. The best fit is A = (13.53 ± 0.09) A, just between 
Fexix at 13.52 A and the Neix intercombination line at 
13.55 A. Second, from a theoretical point of view, we sim- 
ulated the spectra in the region of the Ne ix He-like triplet 
with CHIANTI using the emission measures from table [2] 
and the abundances given in table [3] For the spectral res- 
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olution of the RGS detectors three distinct emission peaks 
are expected. The middle one is always dominated by iron 
lines, even in the case of high densities or strong UV-fields. 
It is always predicted with a photon flux comparable to the 
r line. However, the f line is dominated by emission from 
neon. It is detected with the same flux as the r line in our 
observation (table H] and Fig. [5]), but should be nearly ab- 
sent in the case of high densities or strong UV-fields. From 
these arguments we cannot obtain a numerical f/i value, 
but qualitatively it is clear that the i?-ratio of the Neix 
triplet is high similar to the i?-ratio of the O vn triplet. 

The theory of the He- like triplets is we ll developed 
(|Gabriel fe Jordan! Il969t iPorquet et all l200lh and the R 
and G-ratios are understood. To calculate the expecte d / ji 
we use the formulae given in iBlumenthal et alj l| 1972ft : 

R= L = ^2 m 

i l + ^/^ c + n e /N c 1 ' 

where Rq denotes the low density limit (3.95 for the case 
of O vn). N c and <I> C represent the critical density and pho- 
ton field respectively, where the i?-ratio becomes density or 
radiation sensitive. n e is the electron density. The relevant 
wavelength to excite an electron from the upper level of 
the f to the i line is 1630 A in the O vn triplet and we ob- 
tain the correspond ing stellar flux from an IUE observation 
dValenti et alj|2000ft. We der edden the flux with Ay = 0.25 
( van den Ancker et al.l 11998ft according to the formula of 
iCardelli et alJ (| 1989ft . This in turn gives a radiation tem- 
perature and allows to calc ulate the stellar radiation field 
to get $/<f> c (for details see lNess et al.ll2002D . We use a ge- 
ometric dilution factor W to express the decrease of the 
radiation field with increasing distance r to the stellar sur- 
face i?*: 

Figure[H]shows contours of the expected i?-ratio in the O vn 
He-like triplet depending on the density of the emitting re- 
gion and its position, expressed in stellar radii, where we 
define the stellar surface as r = 1. The region compati- 
ble with the observation is shaded in grey. Qualitatively 
this agrees with the requirement of a large i?-ratio in the 
Ne ix triplet. Clearly we need to look for emission scenarios, 
where the emission region is located above the stellar sur- 
face and the electron number density is below 10 10 cm~ 3 . 

The observation of AB Aur is less well exposed. 
Additionally the i line in the O vn triplet contains two miss- 
ing data bins caused by damaged rows in the detector. The 
f line is clearly present and there can be no contamination 
for oxygen by the near, bright sourc e SU Aur because it is 
strongly absorbed. So we agree with lTelleschi et alj (|2007t ) 
that there are indications for a large f/i ratio, but we can- 
not give a reliable number of counts in the i line. For neon 
we find f/i — 2A^f^ (90% confidence). Due to our small 
extraction region the contribution of SU Aur in the Ne ix 
triplet should be negligible, furthermo re AB Aur seems to 
be iron depleted (jTelleschi et alj|2007ft . but any iron con- 
tamination would further enhance the f/i ratio. This ob- 
servation favours low-density emission regions above the 
stellar surface, too. 



HD 1 63296: VII f/T— ratio 



4 









1 1 1 r. >■•■■! 1 1 r 




1 - 










1 
















: 














_ 






















J 

































7 8 9 10 11 12 

lo 9io( n e) t cm ~ 3 ] 



Fig. 6. O vn He-like triplet f/i ratio depending on the den- 
sity of the emitting region and the position. Shaded are 
those regions of the parameter space within the 90% confi- 
dence range. 
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Fig. 7. Ratio of O vm flux to O vn flux, (colour in elec- 
tronic version only) 

4.5. Cool excess 

iRobrade fc Schmittl (12007D and ICiidel fc Telleschil (|2007ft 
present an analysis of the O vin/O VII ratio as a measure of 
the excess of soft emission in CTTS is respect to stars on the 
MS. These studies show the observed excess to be confined 
to a narrow temperature range about the formation of the 
He-like Ovu triplet at 1-2 MK. In Fig. [7] we add the two 
HAeBe stars with resolved O vn and OviH lines t o thei r 
sample. The MS stars are taken from iNess et al.l ([2004). 
The sample of CTTS has larger total luminosities in oxy- 
gen than the HAeBes. This is a selection bias because the 
analysis requires a clear separation of the O vm and O VII 
line which only grating spectra allow. Those observations 
usually probe the most luminous members of each class. 
The O vm to O vn ratio for a given luminosity is smaller 
for all young stars than for MS stars indicating that the 
younger stars are cooler. Within the CTTS there is a ten- 
tative trend with the youngest stars in the Taurus molecu- 
lar cloud and the Lupus- Auriga star forming region at the 
bottom and older objects like MP Mus and TW Hya closer 
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to the MS. We can plot the ratio only for two HAeBes with 
approximately the same age, but both are set clearly apart 
from the MS. As for CTTS MS stars of the same luminos- 
ity would be considerably harder. In CTTS this soft excess 
is presumably caused by accretion but a contribution from 
winds and outflows is also possible. We present our inter- 
pretation of the cool excess in HAeBes in Sect. 15.41 

5. Discussion 

First, we compare our results for HD 163296 to those ob- 
tained for AB Aur and HD 104237 (Sect. Oil. We present 
reasons to reject the hypothesis that an as yet unde- 
tected companion is responsible for the X-ray emission in 
Sect. HL2J In Sect. 1 5.31 we then discuss the speculation by 



Mq Fe Si S C 



N 



Ne 



ISwartz et al.l (|2005 ) on the basis of imaging Chandra data, 
that HD 163296 may be accretion dominated like the clas- 
sical T Tauri star TW Hya. A full review of many possible 
X-ray generation mechanisms incl uding accretion, winds , 
and disk-related models is given in ITelleschi et "ail (|2007f) . 
Because little information is available to constrain many of 
those models from X-ray observations we refer the reader to 
that publication for a larger variety of models. In Sect. 15.41 
we discuss which mechanism might be responsible for the 
soft X-ray emission and in Sect. 15.51 we present a corona 
from primordial fields as a good candidate for the origin of 
the hard component. 

5.1. Comparison to AB Aur and HD 104237 

HD 163296 and AB Aur share many characteristics. They 
are of the same age and the spectrum differs only by one 
or two subtypes around AO. HD 104237 is of spectral type 
A8. It thus seems natural to expect similar X-ray proper- 
ties from these stars. And, indeed, the total luminosity of 
\ogLx for the first two stars is 29.6 erg s _1 , HD 104237 is 
brighter (logLx ~ 30.5 erg s _1 ), although a contribution 
of about 10% from H D 104237- B has been removed from 
the grating spectrum (jTesta et al.ll2008T) . It is unknown how 
much luminosity the HAeBe primary contributes and which 
fraction is due to the unresolved, close CTTS companion of 
HD 104237. In a fit with two emission components about 
one third of the total emission measure is found at 0.2 keV 
for AB Aur, in contrast to HD 163296 where, according to 
table [H 60 % fall in this temperat ure region. Due to their 
lower signal ITelleschi et al.l (|2007l ) do not split the hotter 
component as we do. Their temperature is comparable to 
our medium component. We find the hardest component in 
HD 163296 to be the one with the lowest emission measure, 
although it is significantly detected in our observation. If 
the same emission mechanism is invoked for AB Aur and 
HD 163296 then a convincing explanation needs to be found 
why HD 163296 is softer. HD 104237 is significantly hot- 
ter, only little emission measure is found below 0.25 keV 
in the reconstructed emission measure distribution. Again, 
the hot plasma might be due to the unresolved CTTS. 

As for elemental abundan ces we left carbon as a free pa- 
rameter instead of nitrogen as lTelleschi et al.1 (|2007l ) did. To 
compare the fitted values for O, Ne, Mg, Si and Fe we first 
need to correct them to the same reference values. On the 
one hand, in AB Aur Ne and Si have the largest abundance 
and O, Mg, and Fe the smallest, in HD 163296 on the other 
hand C, Si, Mg and Fe, those elements which condense on 
grains easily, are found with enhanced abundance. A very 
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Fig. 8. Abundance s for HD 16 3296, AB Au r 
(|Telleschi et al.ll2007l ) and HD 104237 (jTesta et alJl2008t ). 
Error bars show la uncertainties and all abundances are 
normalised to iron. 



similar pattern is found from the spectrum of HD 104237, 
only here S was left as a free parameter of the fit instead 
of C. We show the abundance pattern of all three stars in 
Fig. [S The abundance distribution of AB Aur looks erratic, 
whereas in HD 163296 there is trend with the FIP effect. For 
HD 104237 the pattern is very similar to HD 163296 except 
for a higher neon abundance. However, in young objects as 
HD 163296 with ongoing accretion the physical mechanism 
might be different: The same abundance pattern can also be 
interpreted as an overabundance of elements released from 
refractory grains. The noble gas neon does not condense as 
easily on grains as the other elements in table [3] 

For line emission spectra with only little continuum the 
absolute metallicity is very un certain. Although the lumi- 
nosity's of both stars are equal, ITelleschi et alj (|2007l ) find a 
larger emission measure for AB Aur and at the same time 
a lower metallicity. This likely does not point to a phys- 
ical difference but only shows an ambiguity in the fitting 
process, so we normalised all abundances in Fig. [8] with re- 
spect to Fe. If different emission mechanisms, e.g. a corona 
and a jet contribute to the emission than those abundances 
represent an average value. Therefore we conclude than an 
emission mechanism for soft X-rays, which is possibly based 
in a region of now gaseous former refractory grains, operates 
in HD 163296 much stronger or in addition to the processes 
active in AB Aur and possibly HD 104237. 

5.2. A companion? 

Because HAeBe stars are not expected to posses an outer 
convective layer, X-ray emission detected from these ob- 
jects is often attributed to low-mass companions. The high 
optical luminosity of the main component could overwhelm 
the optical emission of close CTTS. In several systems high 
resolutio n Chandra images allow t o separate multiple com- 
ponents (jStelzer et al.l 2006, 20081) . still as yet undetected 
companions closer to the central HAeBe star might exist. 
In the case of HD 163296 there is no second stellar source 
resolved in the images and the source position matches the 
expecte d coordinates withi n la of the Chandra pointing ac- 
curacy (jSwartz et al. 2005). Furthermore, the same authors 
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obtained HST /STIS data with no sign of binarity. All this 
makes an origin on a late-type companion unlikely. Also, 
the X-ray data differs significantly from both a CTTS in 
the f/i ratio and a weak-lined T Tauri star (WTTS) which 
display acti vity seen in hot plas ma and variable, flaring 
lightcurves (|Preibisch et al.|[2005h . 



5.3. Accretion? 



The r elative softness of the spectrum lead ISwartz et al.l 
(2005!) to speculate that HD 163296 may be dominated by a 
strong accretion component. Accretion of matter could pro- 
ceed along magnetically funnelled streams from the disk to 
the stellar surface. This i s the standard scenario for CTTS 
(jUchida fc ShibatalH98llKoenig]|[l99l . The accretion pro- 
ceeds with free-fall velocity and a strong shock develops on 
the stellar su rface. We fitted HP 1 63296 with an accretion 
shock model (| Giinther et al.ll2007[) and two VAPEC compo- 
nents. In short, our model takes the infall velocities and 
densities as input. It heats the matter up in a strong shock 
and follows the cooling in the post-shock accretion zone in 
a ID geometry. The code explicitly takes into account non- 
equilibrium ionisations. The soft emission around 0.2 keV 
and some of the medium emission could be produced by a 
shock with low densities of 4 x 10 10 cm -3 or less which is 
falling in at around 550 km s _1 . A filling factor / = 7x 10~ 4 
of the stellar surface is sufficient to produce the observed lu- 
minosity. At the given shock parameters about 20 % of the 
total accretion luminosity are emitted in X-rays, the mass 
accretion rate in this case would be M = 10~ n AfQ yr . 
Given the stellar parameters of AB Aur a similar sce- 
nario could apply to that star, too. The accretion rate is 
much lower than the M = 10~ 8 - 1O~ 7 M0 yr " 1 found 
by optical observations (|Garcia Lopez et al.ll2006l) or mod- 
elling of the spectral energy distribution fSE DI dDent et al 



2006). This resembles our resul ts for CTTS (|Giinther et al 



2007; lGimther_&_Schmitt 2007) and we attribute it to non- 



uniform accretion spots, which only produce X-rays in the 
innermost hottest and densest parts o f the spot. Our ar- 
gument contrasts iTelleschi et al.l (|2007f) : They take the op- 
tically determined accretion rate and the density and dis- 
tribute it over the stellar surface, which results in a flux 
much larger than observed. Our models are more sophis- 
ticated and explicitly give the fraction of the accretion lu- 
minosity that is radiated in the soft X-ray band, but they 
are invalid in the case of strong absorption. The low column 
density found in the observation supports the validity of our 
model. We expect the post-shock cooling zone to penetrate 
the stars down to 0.97i?*, but most of the X-ray emission is 
generated close to the stellar surface. Due to the low-density 
of the pre-shock material the optical depth is low and the 
X-rays escape. Observationally a large optical depth would 
be seen as a significantly suppressed resonance line in the 
He-like triplets, whereas the intercombination and the for- 
bidden lines remain unaltered. In all stars analysed so far 
the strength of resonance line is compatible with the ex- 
pected G-ratio, thus proving that the emission is optically 
thin. 

In addition to the energetics the accretion shock model 
explains the small f/i ratios found in CTT S with the high 
densi ties in the post-shock cooling zone (| Giinther et al.l 
2007), but in HD 163296 and AB Aur we find an f/i which 
is fully compatible with the low-density limit. This excludes 



regions close to the stellar surface as emission origin because 
here the stellar radiation field of an A type star would radia- 
tively shift emission from the f to the i line, if the accretion 
zone is not shielded from the radiation. If it is, e.g. by the 
outer layer of the accretion stream, then the X-ray emission 
should be more absorbed, but the measured absorbing col- 
umn density agrees well with the stellar optical reddening 
fSect. FO]) already without an ext ra absorption compo nent. 
So we reject the original idea by ISwartz et al.l (|2005| ) that 
the soft X-ray emission in HD 163296 is powered by accre- 
tion shocks. 



5.4. The origin of the soft component 

The soft radiation in our three component model originates 
in a plasma with temperatures around 0.21 keVw 2.4 x 
10 6 K. Plasma can be heated to this temperature if gas 
moves with and passes through a strong shock (it can be 
slower then the 550 km s^ 1 mentioned above because we 
try to explain the soft component only), according to the 
Rankine-Hugoniot conditions, which lead to the following 
formula: 



^shock 



400 km s" 



kT 



0.21 keV' 



(3) 



where we transform the kinetic energy of the pre-shock ve- 
locity Vshock in the shock rest frame into the post-shock 
thermal energy kT. Other situations than the accretion 
shock discussed above can give rise to such scenarios. 
Shocks occurring in winds or outflows easily fulfil the condi- 
tions on the density and distance to the photosphere set by 
the analysis of the He-like tiplets. In AB Aur temporal vari- 
ation in the X-ray luminosity was detected which is compat- 
ible with the period of 42 h, observed in Mg II and He I lines 
which are tracers of winds and outflows (Tell eschi et al.l 
120071 ) . Shocks could develop where fast and slow winds col- 
lide either through magnetic collimation or because of tem- 
poral variability of the launching velocity. Here we wan t to 
discuss a scenario not covered by ITelleschi et al.l (|2007f ) . 

HD 163296 drives two powerful collimated jets which 
are visible in corqnographic images obtained with HST 
(|Devine et al.ll2000l : lGradv et al.ll2000D . In this respect it is 
similar to the lower mass CTTS DG Tau , where also X-ray 
emiss ion has been detected from the jets (jGudel et al.| [2005. 
2008) . The high extinction of the central source DG Tau it- 
self allows to detect a signifi cant offset between the sof t 
and the hard X-ray emission ([Schneider fc Schmittl [2008) , 
because it absorbs all contributions from the star in the 
low energy band. This offset points to an origin of the soft 
emission in the jet forming region, possibly as a collima- 
tion shock or in an emerging kno t where faster electa catch 
up with previous older outflows (|Giinther et al.ll2008f ). We 
searched for an offset in the Chandra data for HD 163296, 
but could not find any deviation from the expected PSF, 
so the emission region must be within 100 AU of the star. 
Line ratios of Fen indicate electron densities between 10 4 
and 10° cm~ 3 in the knots of the jet (|Wassell et al. 2006). 
Presumably the density increases towards the star, so we as- 
sume a pre-shock particle number density of uq — 10 6 cm -3 
in the shock region. With this assumption we can apply 
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the approximation to the post-shock cooling length d coo i of 
iRaga et al.l (|2002[ ): 



d coo i « 0.8AU 



10 6 cm~ 3 
n 



^shock 



400 km s" 



4.5 



(4) 



and find that d C ooi is only 0.8 AU. Dividing the volume 
emission measure (VEM) by the density we obtain the vol- 
ume and - using the equation above - we can transform 
thi s to a radius R of th e shock front (in AU) (for details 
see iGiinther et aTll2008D 



observations bv lHubrig et al.1 (|2007[ ) support this idea with 
signatures of magnetic fields in Can H and K lines likely 
of circumstellar origin. Unfortunately the signal does not 
allow a qualitative measurement of the field strength which 
is needed to estimate the radius of the magnetosphere for 
alternative model s of magnetically confined winds. Also 
Wa de et al.l (|2007f ) have questioned the method used in the- 
ses studies. 



10 b cnT 



0.5 



VEM 
10 52 cm -3 



05 '0.21keV\ 1125 



kT J 



where we assume a cylindrical geometry with the shock at 
the cylinder base. The mass flux Mshock through the shock 
is independent of the density, because smaller densities lead 
according to Eqn. \5\to larger shock areas: 



hock 



6 • 10 



_ U M / VEM \/0.21keV 



V io 52 



cm 



kT 



1.75 



(6) 



For HD 163296 a mass outflow of the order M s hock ~ 
10~ 10 Mq yr _1 is sufficient to power the soft X-ray com- 
ponent, this is two orders of magnitude smaller than the 
mass loss rates determined from UV and optical data 
(|Deleuil et alj|2005t IWassell et al.ll2006fh Only a small com- 
ponent of the total jet needs to be fast enough to be shocked 
to X-ray emitting temperatures. 

This scenario seems feasible in the light of the opti- 
cal observations. The velocity of the jet determined from 
a combination of radial velocity and proper motion is 



360 ± 82 km s" 



yelocity and proper motion is 
(jWassell et al.l I2006D where the 



radial velocity is taken as the average of a broader distri- 
bution at the position of the innermost knot. The offset 
between the central source and the knot A is 9", at the dis- 
tance of HD 163296 this corresponds to 1100 AU. Certainly 
the jet is launched faster and decelerates as it goes through 
the shock and later moves through the surrounding matter. 
In the shock rest frame the velocity v\ drops to a quarter of 
the pre-shock velocity vq, where the matter passes through 
the shock front: 4 • Vi = Vq. To heat the plasma to the re- 
quired temperatures we need Av = vo — V\ = 400 km s . 
This is all fullfilled if the shock front travels outward along 
the jet at Wf ron t = 230 km s _1 in the jet rest frame. The jet 
emerges at t>i aU nch = Wf ron t + i>o = 750 km s _1 and the knot 
moves at Wj e t = Vfront + v% = 360 km s _1 with respect to 
the star. 

Typical jets o f CTT S reach up to 400 km s" 1 
(|Eisloffel fc Mundtl 11998), co mparable to the case of 
HD 163296. iRodriguez! (jl995T) report velocities of up to 



5.5. The origin of the hard component 

To heat the plasma emitting the hot component shocks with 
(5) velocities above 1300 km s _1 are needed but the optical ob- 
servations show no lines which are shifted by this amount. 
Other heating mechanism need to be responsible and a 
promising candidate are reconnection events in magnetic 
fields in analogy to the solar and other late- type coronae. 
No polarisa tion is found in photospheric lines (|Hubrig et al.l 
2006, l2007t ) but this may be due to strong small-scale fields 
with only little global dipole contribution. Young A stars 
may either still posses primordial magnetic fields or gen- 
erate t here own. In this evolutionary stage the simula- 
tions of ISiess et al.l (|2000f ) predict HAeBes to have a thin 
outer convective layer where a dynamo could operate or the 
shear in rapidly rotat ing stars could cause magnetic fields 
(|Tout &: Pringld[l995D . Given the detection of polarisation 
in circumstellar matter and the collimation of the outflow 
we expect that magnetic fields exist on HD 163296 itself. 

No flaring activity has been unambiguously observed 
on a HAeBe star and especially AB Aur, HD 104237 and 
HD 163296 show only little varia bility in their hard com po- 
nents during the observation, but lGiardino et all i|2004h can 
at least give good indications that the flare they observed 
occurs on a HAeBe star and not on the known, but unre- 
solved CTTS companion. We postulate that the medium 
and hard component of the X-ray emission originate in a 
stellar corona. At high temperatures the contribution to the 
O VII and the Ne ix He-like triplet is small so the observed 
f/i ratio cannot trace the UV field in the emission region 
of the hard photons. 



1000 km s out in the jet, so an initial velocity of 
750 km s _1 is not unreasonable. This X-ray generating 
mechanism may operate in addition to any of those already 
suggested for AB Aur or HD 104237. Possibly HD 163296 
is softer because in AB Aur has no jet which could emit 
soft X-rays. 

Several models for the acceleration and collimation of 
outflows have been suggested even for the class of CTTS 
where the sample size is much larger than for HAeBes, still 
the exact collimation mechanism is uncertain. Yet, all mod- 
els rely on a magnetic field. So the existence of two jets 
already makes the presence of magnetic fields in the en- 
vironment of HD 163296 very likely. Spectropolarimctric 



6. Conclusion 

We present the observation of HD 163296 with XMM- 
Newton. This is after AB Aur the second grating spectrum 
of a HAeBe, without -as far as we know- an unresolved 
companion. The signal-to-noise ratio (SNR) is significantly 
better than in the first observation. Congruently both stars 
show soft spectra with little temporal variability. In all ana- 
lysed He-like triplets the line ratios point to an emission ori- 
gin in regions of low densities and weak UV photon field, 
that is at a few stellar radii. On HD 163296 we find a FIP 
abundance pattern typical for inactive stars. 

Both the spectral shape and data from other instru- 
ments make an origin of the X-ray in an unknown com- 
panion to HS 163296 very unlikely. Therefore we believe to 
present genuine X-ray data for HAeBe which can be repre- 
sentative for the whole class of this objects for all we know 
today. 

The optically observed mass flux rates are much larger 
than required by the energetics of the X-ray radiation in 
the case of accretion shocks or emission from outflows, but 
the large //i-ratios rule out an accretion shock scenario 
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on the stellar surface. We favour a model where the main 
contribution to the soft X-ray emission comes from shock 
heated plasma in the jet or wind. Possible geometries in- 
clude shocks in unsteady winds or collimation shocks on the 
base of the jet. The collimation of the jets is likely related to 
a magnetic field, whose origin is unclear. This field can be 
a frozen-in primordial magnetic field or the star generates 
one in a thin convective layer. We find no other convincing 
explanation for the hard emission than heating by ener- 
getic events from the magnetic field. Thus we predict that 
with longer observation times we should observe flare-like 
behaviour on HAeBes. 
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50OR0105. 



Appendix A: Error estimation for line ratios 

For experiments with low count numbers the observed 
counts can be taken to be Poisson-distributed. For longer 
exposure times or stronger sources the Poisson distribution 
can be reasonably approximated by a Gaussian distribu- 
tion. In this case the error on a line ratio can be calculated 
according to the conventional error propagation formulas. 
For weak sources, however, the Poisson distribution leads 
to significantly asymmetric errors, which are dominated by 
the statistical uncertainty in the weakest line. 

We set up a Monte-Carlo simulation to determine con- 
fidence intervals. The instrumental and source background 
is combined in a total background b, which we assume to 
be uniform over a range of a few A. It is measured from 
a larger line-free region on both sides of the triplet of in- 
terest. Thus the mean background level is well known. In 
the case of the O vn triplet in HD 163296 we determine the 
background from the line-free regions in the range 20-25 A 
to b = 30 counts A -1 . 

For a fit with i lines, the spectrum is divided into i bins 
centred on the lines, each bin is AA^ wide, thus the expected 
background per bin is 

Bi = b-AXi (A.l) 

The line profile is given by the instrumental profile, in 
the case of XMM-Newton the data is best fit by a Lorentz- 
profile. For each bin i we compute the fraction of the 
total line spread function of line j that is expected in bin 
i. For known line fluxes fj then the observed number of 
counts in each bin can be calculated. We write this in form 
of a matrix equation 

L = A F + B, (A.2) 

where A is constructed from the atj, B from and F 
from fj. L represents the observed counts per spectral 
region. Since the sum of Poisson distributions is again 
Poisson-distributed, we can take the number of counts in 
each bin as Poisson-distributed, even it is calculated as 
the sum over several detector bins. Thus the values of 
L are Poisson-distributed and we simulated 100 000 real- 
isations of a Monte-Carlo calculation for each bin, taking 
the observed count number as expectation value. Inverting 
Eqn. IA.2l gives 

F = A- 1 -{L-B). (A.3) 




Fig. A.l. Statistical distribution of the i?-ratio based on 
10 5 realisations of a Poisson statistic for the forbidden and 
intercombination line. The red/grey line shows the cumu- 
lative distribution function. 



From the simulations we use this to obtain the values of the 
individual line fluxes /. Negative numbers occur, if the sim- 
ulated value of the Poisson distribution happens to be below 
the fixed background value. We reset these numbers to 0; 
in the analyses of real data they represent situations, where 
the line seems to be absent, thus producing extreme values 
for the line ratio. The simulated count ratios are scaled by 
the effective area and the flux ratios are calculated. The 
simulated distribution of the //i-ratio is shown in Fig. IA.ll 
From the cumulative distribution (red/grey line) the confi- 
dence limits can be read off. 
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Fig. 3. High resolution RGS spectra with the best-fit model from the global fit in Sect. I4.2I with the corresponding fit 
residuals. The data is binned to a minimum of 5 counts per bin. black: RGS1 red/grey: RGS2 



